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We study the distribution of the phase angle and the magnitude of the fermion determinant as 
well as its correlation with the chiral condensate and the baryon number for QCD at non-zero 
quark chemical potential. Results are derived to one-loop order in Chiral Perturbation Theory 
(ChPT), as well as by analytical and numerical calculations in QCD in one Euclidean dimension. 
We find a qualitative change of the distribution of the phase of the fermion determinant when the 
quark mass enters the spectrum of the Dirac operator: it changes from a periodicized Gaussian 
distribution to a periodicized Lorentzian distribution. We also explore the possibility that some 
observables remain weakly correlated with the phase of the fermion determinant even though the 
sign problem is severe. We discuss the practical implications of our findings on lattice simulations 
of QCD at non-zero baryon chemical potential. 
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1. Introduction 

It is well known that standard lattice simulations based on importance sampling are not pos- 
sible for QCD at non-zero chemical potential, since the measure in the functional integral is no 
longer positive definite. In the past years a few methods have been proposed to circumvent this 
problem, at least for small values of \ijT [1]. They all share a common philosophy: choose wisely 
a simulation ensemble - either at zero or non-zero imaginary chemical potential and fluctuations 
will allow the exploration of the target ensemble at non-zero baryon density. In a natural way, we 
are thus led to consider the overlap between the simulation ensemble and the target ensemble, and 
to devise strategies to maximize it. This note, which is based on a more extended publication [2], 
reports our first results towards this goal. 

Fluctuations increase with temperature so in general ensembles generated a higher tempera- 
ture allow a better re- weighting [3], an easier calculation of the Taylor coefficients [4], or a safer 
extrapolation from imaginary chemical potential [5, 6, 7, 8]. Further, the shape of the distribution 
function plays a major role in the density of states method [9, 10] and some assumptions on its 
shape are needed. 

The main goal of our investigation is to put these heuristic considerations on quantitative 
grounds by a combination of analytic studies and numerical simulations of chiral perturbation 
theory and QCD in one dimension. 

We will base our discussion mostly on the 6 -distribution (8(6 — 6'))^ f and the constrained dis- 
tribution (0 8(6 — 6')} (the phase of the fermion determinant is denoted by 6'). The 6 -distribution 
will help assessing the overlap between the simulation and target ensembles, and the constrained 
distribution shows how averages are built up in the spirit of the density of states method (the inte- 
gral over 6 obviously gives the full expectation value (&)). The distribution of an observable over 
the phase allows us to address which range of the phase is essential for the full expectation value 
of G. 



2. Results from Chiral Perturbation Theory 

Let us first recall that the 6 distribution in the full theory 

V V " 1 J dA\det(D + H7b+m)\ Nf e lN f e e- s ™ 

has a simple expression in terms of the phase quenched distribution 

(8(6 - 6')) Nf = e ieN f^(8(6 - 0')) m . (2.2) 

^N f 

In the following, we will mostly consider the two flavor case (Z\ + \* = Zm and Zi+i = Z2): 

(8(6 - 0')>i+l = e mZ ^(8(6 - 0')>i+i*. (2.3) 
^1+1 

The calculations follow the general method of [11]: the 6 distribution is obtained from the 
moments of the phase factor 

(S{0-0')) N , = — £ e-^(e'P e ') Nf . (2.4) 

Z71 p=-oo 
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Figure 1: The real part of the distribution of the phase (8(9 — 0'))i+i (solid curve) for a small AGq = 0.2 
(left) and a large AGq = 10 value (right). Also shown is the real part of the distribution of the baryon number 
over 9 (dashed curve). In either cases the 9 distribution is normalized to 1, while the integral of distribution 
of the baryon over 9 is zero. This directly illustrates the severity of the sign problem at large Go 



First we consider /J. < m n /2 where the leading order difference between the phase quenched 
and full free energies is determined by the one-loop corrections. The 6 distribution to this order is 
given by a periodicized Gaussian 

(8(6 - 00)1 + 1 = V e -( e + 2 ™) 2 /AGo+AG _ (2.5) 

where Go is the free energy difference between neutral and charged pions. The distribution of the 
baryon number ub over 8 is given by 

(n B 8(6-e')) l+l = (lim^AG (-n,fl)) £ (1 + i d ±^ ) ^^ e -(e + 2^/AO 0+ AG _ 

(2.6) 

Note that, since (ns) = in ChPT, the integral of the above distribution is zero as it should In 
contrast, the integral of 6 -distribution is unity. The similarity of the two distributions displayed in 
Fig. 1 clearly shows that the details of the distribution are important to get the correct results for 
ns, especially when Go is large, i.e. when the sign problem is severe. 

In order to make contact between these calculations and lattice studies, we have also have 
computed the phase quenched as well as the partially quenched 6 distribution where the ensemble 
is generated at zero chemical potential [12]. The relationship between the 6 distributions derived 
in one-loop ChPT is surprisingly simple: Whether we compute the width of the Gaussian for the 
6 -distribution in the full ensemble generated at fi, or the partially quenched ensemble generated at 
H = 0, or in the quenched ensemble, or the phase quenched ensemble, we find exactly the same 
result. These results directly apply to the interpretation of lattice studies at small chemical potential 
[12, 13]. 
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We now proceed to the computation of the -distribution for 11 > m n /2. Only the even mo- 
ments are accessible in ChPT and give the distribution 



(5(20-20')}^ = 1 £ e-^V'^V (2.7) 
The leading order contribution now enters already at mean field level. This leads to the result 

(5(28-28')),+, = ^ f*> (2.8) 

7T COSh^VLfi) — cos(20) 

where L# is the difference of the free energy densities in the full and the phase quenched theory. 
The quenched result is obtained simply removing the general factor e 2l6 e VLB . In both cases the 
distribution is a periodicized Lorentzian. This is in sharp contrast with Gaussian obtained at low 11, 
and implies that the hypothesis leading to the central limit theorem are not realized for /I > m^/2. 
Further, we discuss the the distribution of 

F = |det(D + i U7o + m)|/det(D + m) = exp(/). (2.9) 

This quantity was studied in lattice QCD [12] using the Taylor expansion method. We have verified 
the assumption [12] that the distribution remains Gaussian even for a fixed value of F. Of course 
this applies only to jU < m K /2. 

Since determinants fluctuate by many orders of magnitude we feel that it is more appropriate 
to analyze the distribution of logF = / instead. The two distributions are related by a simple 
transformation 

(8(f-f')) = F(S(F-F')). (2.10) 

The fluctuations of / are induced by the gauge field fluctuations. 

We find that for ll < m n /2 to one-loop order in ChPT the distribution of / is Gaussian 

_ (f-N f E f /2) 2 

W-f'))N f = -^=e 2 °f . (2.11) 

where we have given the result for an arbitrary number of flavors Nf. Note that Ef and ay do not 
depend on Nf, see [2]. 



3. QCD in one Euclidean dimension 

Consider now one-dimensional QCD with a massive staggered fermion and gauge group 

U(N C ) 

Z Nf (Hc,H)= / dUdetM (3.1) 

JU(N C ) 

where 

detM = 2-" Nc det[e nlic +e- n ^ + e n ^U + e- n ^U f ], U € U(N C ) (3.2) 

and dU is the Haar measure. The group integral can be solved analytically [14, 15], which shows 
explicitly that the partition function does not depend on fx, as expected since there are no baryons 



4 



Fluctuations, correlations and the sign problem in QCD 



M. P. Lombardo 






P =i - 

p = 2 
p = 3 
p = 4 




p = 5 . 







Figure 2: 1/(3) lattice QCD in one dimension for m q = 0.1 i.e. jj. c = sinh -1 m q . Left we show the logarithm 
of the p— th moment as a function of p, for = (0. 12, 0.28) > /i c , showing the expected linear behavior Eq. 
(3.4). The moments themselves as a function of jj. for p = 1 (the smoothest) to p = 5 (right panel), show a 
sharp change at jti = fi c leading to the qualitative modifications of the 9 distribution function described in the 
text. Note also the different behavior of even (black symbols) and odd (red) moments, respectively (right). 



in a U (N c ) model. However, the quenched model has a transition to a pion condensed phase at 
/J. = /l c . So it resembles ordinary QCD with p, < jIIb/3, in contrast to SU(N C ) in one dimension 
where the transition of the quenched and the full model are very close to each other (which is of 
course the reason why we decided to work in U (3) rather than in SU (3)). 
Again the 6 distribution follows form the moments 



(e 2 'P e = / dU——. (3.3) 
x 1 Ju(n c ) faxPM 



For p <»ij/2 we obtain 



,2\P 2 



(&*) = (^*)?=(\-^\ . (3.4) 
With Q. = — log(l — p 2 /p 2 ) this yields the quenched distribution 

(8(6-d')) = ^=e- e2/a for ju < He, N c -»• oo, (3.5) 
v 7r£l 

valid for 26 G [—00,00]. 

For p > /J, c the even moments can be computed exactly for any N c . They are given by 

= e - 7 "WcV i (3 . 6) 

and the quenched 6 distribution again becomes a periodicized Lorentzian 

(8(26 -2d')} = - ^ inh (2^) for M>/Ic 2 0€[-W,4 (3.7) 

x v ;/ 7rcosh(2«A^ ci u)-cos(20) ^ ^ ' 

The results have been verified numerically with a good accuracy, see Figure 2. 
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4. Summary, and comments 

We have studied the distribution and related properties within chiral perturbation theory 
and one dimensional QCD. The qualitative features of the results are the same. In either case 
the baryons are omitted by fiat, and despite the absence of baryons there is a highly nontrivial /I 
dependence of these distributions. 

Our main observation concerns the non trivial changes of the distribution when the chemical 
potential /J. exceeds the quenched threshold for pion condensation pL c = m K /2, it can be summarized 
as follows: 

• {J. < m n /2 , i.e. quark mass outside the eigenvalues: Gaussian distribution. 

(8(6 - d')) = —== V e -(e+2^) 2 /AG _ 

• /I > m n /2 , i.e. quark mass inside the eigenvalues: Lorentzian distribution. 

MO 7() '\\ - 1 sinh(VL B ) 
(5(20-20 )) - - cosh(VLB) _ cos(20) - 

We reiterate that jl c = m n /2 is not a critical point of the full theory [16], and any shaip 
change at this unphysical threshold implies almost automatically specific difficulties in numeri- 
cal approaches relying on extrapolation from low chemical potential values. 

Further comments concern the validity of the central limit theorem: for jU < m K /2 the distribu- 
tion is Gaussian, thus fulfilling the conditions of the central limit theorem, and in agreement with 
the behavior found in lattice simulations. However, the Lorentzian shape of the distribution of the 
phase valid for larger values of the chemical potential shows that one should not take for granted 
that the conditions for the central limit theorem are satisfied. 

The analytical results show that exponentially large cancellations may take place when in- 
tegrating over 0, needed to measure correctly the baryon number. The extreme tail contributes 
significantly to the results. A small non Gaussian term in the tail of the -distribution therefore 
could be the dominant term after integration over 0. The precise form of this tail is of course 
difficult to access numerically. 

In QCD in one Euclidean dimension the same behavior has been observed by a direct evalua- 
tion of the involved partition functions, either analytically or numerically. We have studied in detail 
the behavior of the moments of the distribution. The change of the distribution from Gaussian to 
Lorentzian is clearly seen in the sharp change of the moments at /i c . In the same simple model, the 
distribution of the gluonic observables can be studied numerically. The results for the distributions 
of the plaquette - to be reported elsewhere [20] - can be used to demonstrate the effectiveness of 
the Fodor-Katz re-weighting [3] : we have explicitly shown that the use of configurations generated 
at high temperatures maximizes the overlap with a cold, dense target ensemble. 

We close with a few general comments: the validity of the Gaussian assumption of the distri- 
butions are limited }l < m K /2: Hence, the point /I = m n /2 looks more dangerous after this study: 
qualitative changes of the distribution, discontinuity of moments. 
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This rather pessimistic outcome might be partially mitigated by two observations: First, the 
6 distribution might be particularly difficult - experiments with plaquette or other observables in 
simpler models might be useful and give a different insight. Second, and perhaps more significantly, 
the results of this study do not include baryons. The physical motivations for the success of any 
practical method for circumventing the sign problem relies on baryonic fluctuations [5, 6]. One 
way to see this is to consider the average phase factor from Taylor expansion [4] : 

(e 2ie } 1+1 * = e L3T ^-^ 2 (4.1) 

or, perhaps more simply, its log-derivative: 

^-log(e 2,e ) 1+1 * = ^l°gZi + i - ^-logZi+i. oc (n B (n)-ni(n)) (4.2) 

which equals the off diagonal susceptibility, and indicates an interplay between the sign problem, 
the temperature, and the baryonic contributions [17], [18]. Current and future work will address 
this issue [19, 20]. 
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